Isospin breaking of the narrow charmonium state of Belle at 3872 MeV as a deuson 
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The narrow charmonium state near 3872 MeV reported by the Belle collaboration and confirmed 
by CDF lies almost exactly at the D°D*° threshold. As was predicted many years ago it can be 
a DD* deuteronlike meson-meson state called a deuson. If so, its spin-parity should be 1 ++ , or 
possibily , and its isospin predominantly 0. An important decay mode should be via D°D*° to 
D°D tt . Large isospin breaking is expected because of the isospin mass splitting between the neutral 
and charged D (D*) mesons. Because of this large isospin breaking the decay X(3872)^ J/^pP would 
be allowed, while J/ipa would be forbidden by C-parity, as indicated by the present data. 

"■ PACS : 12.39.Pu; 14.40Aq; 14.40Gx; 14.40.Lb 

Recently the Belle collaboration [1] reported a new narrow charmonium state at 3872 ± 0.6(stat) ± Q.5(syst) MeV 
£NJ , and with a width < 2.3 MeV (95% C.L.), which is consistent with their resolution. This new state, X(3782), has been 
confirmed by the CDFII collaboration at Fermilab [2] (with a mass 3871. 3±0.7±0.4 MeV). It is 60-100 MeV above the 
expected spin 2 cc ( 3 -D C 2) state [3,4]. Belle sees the new state in the ir + n~J/ip invariant mass distribution from B^ 
j>! ■ decays produced in e + e~ collisions at the T(4) resonance. They find 34.1±6.9±4.1 events of B — > K^tt+tt^ J/tp and a 
10.3fi signal significance for the observed resonance peak. CDF has 730±90 X(3872) candidates in the J/^-k + tt^ mass 
on : distribution produced from initial pp collisions. Many recent papers have discussed possible quantum numbers and 
interpretations of this new state [5-15] finding that interpreting the X(3872) as a conventional cc state has difficulties. 
, On the other hand, a multiquark or 2-meson composition, whose existence have been suggested long ago [16-18] seems 
more likely. The fact that it is very narrow, excludes many states which should decay strongly to DD. Furthermore 
f* — \ the X(3872) mass is is not near any expected charmonium state and the nonobservation of X — > Xcil suggests C = +. 
C*~) ■ Therefore the X(3872) may well be a non-gg state. 

Here we discuss the possibility that X(3872) is a DD* deuteron-like state or a deuson. After Belle announced their 
first note on their observation on the X(3872) [1] I wrote a short reminder [19], for the web, that such a state had 
been predicted long ago with specific quantum numbers near the DD* threshold together with many others near the 
D*D* thresholds and about 50 MeV below the BB* and BB* thresholds. See Table 1. Here we elaborate on the 
possibility that X(3872) is such a deuson and study in particular the large isospin breaking for such a state. 

No doubt the deuteron is a multiquark state which to a good approximation can be understood as a proton-neutron 
system bound by mainly pion exchange. It is therefore natural to ask [20] the question: For which quantum numbers 
O i| of two-meson composites is the well known pion exchange mechanism attractive and comparatively large? This was 
found to be the case for several light meson-meson channels with quantum numbers where problematic resonances 
have been seen, but the attraction was found to be not strong enough to alone bind two light mesons such as KK* . 
Then extending the model to the heavy meson sector [21] using a similar simple nonrelativistic framework as for the 
deuteron [22] many deuteronlike states were found to be bound or nearly bound. These are listed in table 1. At 
that time similar arguments were presented by Manohar and Wise [23] and Ericson and Karl [24] , although not with 
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■ explicit model calculations. 
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1++ 


w 3870 
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1++ 


« 10562 
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0++ 


w 4015 
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0++ 


w 10582 
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« 4015 


B*B* 


0"+ 


w 10590 


D*D* 


1+- 


« 4015 


B*B* 


1+- 


w 10608 


D*D* 


2++ 


« 4015 


B*B* 


2++ 


w 10602 



TABLE I. Predicted masses [21] of heavy deuteronlike states called deusons. These are close to the DD* and the D* D* 
thresholds, and about 50 MeV below the BB* and B* B* thresholds. All states have 1=0. The mass values were obtained from 
(a rather conservative) one-pion exchange contribution only. 
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The X(3872) looks very much like one of the two first deuteronlike DD* states predicted near 3870 MeV (see table 
1 which is from table 8 of Ref. [21]). The states are of course all eigenstates of C-parity, or (\DD* > ±\DD* >)/\/2 
although we here denote them DD* (and in the figures even without the "bar" DD*). We note in particular: 

• Its spin-parity should be 1 ++ or possibily h . For other quantum numbers pion exchange is repulsive or so 
weak that bound states should not be expected. 

• No DD nor BB deusons are expected since the three pseudoscalar coupling (in this case the DDn coupling) 
vanishes because of parity. 

• If isospin were exact the Belle state as a deuson would be a pure isosinglet with a mass very close to the DD* 
threshold. For isovector states pion exchange is generally one third weaker than for isoscalar states. Therefore 
all predicted states were isosinglcts. But as we shall see the isospin breaking will be large. 

• The observed peak (at 3872 MeV) is within error bars at the D a D*° threshold (3871.2 MeV), while the D+D*~ 
channel, 8.1 MeV higher, is closed by phase space. 

• As a deuteronlike state with small binding energy (for the deuteron it is 2.22 MeV) the Belle state should be 
large in spatial size. It should then have a very narrow width since annihilation of the loosely bound DD* state 
to other hadrons is expected to be small, although states containing the J/ip are favoured compared to states 
with only light hadrons due to the OZI rule. 

• A large part of its width should be given by the instability of its components or the D* widths, corrected for 
phase space effects. An important decay should be via D a D*° component to D°D°Tr a since the other charge 
modes lie about 2 MeV above the resonance. Using isospin for the decay couplings and the width measurements 
[25] for D* + together with the dominance of the D°D*° component estimated below, one finds that the width 
to D°D°tt° should be of the order 30-40 keV, and the D° D° 1 / D° 'L> 7r° ratio 62% . 

The deuteron 

Let us now discuss our model for a deuteronlike X(3872) in little more detail. We begin by recalling the essential 
facts for understanding the deuteron. There pion exchange is the dominant mechanism which determines the long 
range part of the proton-neutron wave function and most of the binding energy. It is important to remember that, as 
for the deuteron, the potential comes from an axial dipole-dipole interaction where the tensor potential with its 1 jr 
term (T(r) below) is very large (although it must be regularized at short distances). Without the tensor potential and 
the associated D-wave the deuteron would not be bound. The cental potential C(r) provides only about one third of 
the binding energy, while the freedom for the deuteron to flip the nucleon spins and make transitions between a 3 5i 
and a 3 Di state through the tensor force lowers the energy enough to be bound. This is important in spite of the fact 
that the D-wave interaction is repulsive, and the fact that the D-wave has an angular momentum barrier. 

The central and tensor parts of the potential are given by the functions 

C(r) = , (1) 

T(r)=C(r)[l + — + -r-^] , (2) 

As can be seen from these expressions for distances below a few fermi the tensor potential T(r) is an order of magnitude 
larger than the central C(r), and can therefore not be neglected. In matrix form the deuteron potential can be writtem 
compactly in the form 



- -fvb 



(3) 



where the first matrix index refers to the 3 Si component and the second index to D\. The factor Vq = j^f^ w 1.3 

MeV is fixed by the pion- nucleon coupling constant (g), the pion decay constant, and by the constant 25/9 which 
comes from the spin-isospin factor for the 3 Si state. 

Another fact one should remember is that one cannot predict exactly the deuteron binding energy, without a very 
good detailed knowledge of the very short range interactions. The actual measured value, 2.22457 MeV, is a net effect 
of different contributions of opposite signs, most importantly from the repulsive kinetic term and the attraction from 
all the potential terms. One cannot predict its value by scaling arguments from the other scales in the problem. But, 
given the pion- nucleon coupling constant, and a regularization model for the tensor potential, one can say that two 
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nucleons is most likely to be bound, mainly by pion exchange, and with a small binging energy for 1=0 and for a 
specific mixture of 3 Si and 3 D\. For other quantum numbers the binding is much weaker or repulsive, and no bound 
states exist, (although the 1 5o state is nearly bound). 

But once the binding energy is known and used as an input, to fix the regularization of the tensor potential at short 
distances, the dominant long distance parts of the wave function of the deuteron can be predicted. See Fig. 1. We 
use the same regularisation as in [21] with a parameter A, which can be related to a spherical pion source around the 
hadron with rms radius %/lO/A = 0.624fm/[A/GeV]. From the wave functions one can then predict static properties 
of the deuteron such as the quadrupole moment in agreement with experiment. 
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FIG. 1. The deuteron wave functions 3 Si and 3 Di, when the the pion exchange potential is regularized at small distances 
to get the right binding energy. Apart for very short distances this agrees with more detailed models. 

The X(3872) as a DD* deuson 

This same argument as for the deuteron hold of course true also for deusons. One can determine for which quantum 
numbers pion exchange is attractive and strong enough so that bound states are expected, but hardly the exact 
binding energy. But from the strength and sign of the potential one can exclude already a large number of possible 
deuson candidates with given flavour, C-parity and spin-parity quantum numbers. In fact only for J PC — 1 ++ and 
possibily h and for 1=0 (in the exact isospin limit) can one expect such states to exist just below the DD* threshold, 
as was found in [21]. 

For a 1 ++ DD* system one has a rather similar potential as for the deuteron with an 3 Si and a 3 Di-wave. It can 
be written in a similar form as Eq.(3). One finds in the exact isospin limit a potential 




regularization parameter A is here 1.355 GeV and the isospin splittings between charged and neutral D, D* and n are neglected, 
by using the average D, D* , tt masses in the calculation.) 
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where the overall constant 7 = 3 is determined by the spin-isospin and C-parity of the state [21]. One finds assuming 
a 0.5 MeV binding energy the wave functions shown in Fig. 2. 
Isospin breaking 

However, since the binding energy of a such a deuteronlike state is expected to be at most a few MeV or of the same 
order as the isospin mass splittings one must expect substantial isospin breaking in the wave functions. For a pure 
1 = state one has equal contribution of the two components D°D*° and D+D*~ in (|D L>* > + \D+ D*~ >)/\/2, 
but since D°D*° is 8.1 MeV lighter than D + D*~ it should have a greater weight than D + D*~ . In other words this 
means that there is a substantial 1=1 component in the state. Now, assuming exact isospin for the pion couplings, 
the \D°D*° > and the \D + D*~ > channels are coupled through the matrix 

(11) 09 

The eigenvalues of this matrix are 3 and -1. The first corresponds to 7 = 3 in Eq.(4) i.e. to the 1=0 linear combination, 
where one has attraction, while the second (-1) corresponds to the weaker and repulsive 1=1 DD* channel. 

Now to estimate the isospin breaking we take account the fact that the charge modes of the D , D* and the pion 
are split in mass. Therefore (a) the kinetic terms for \D°D*° > and \D + D*~ > are slightly different, (6) the potential 
is slightly different for tt° exchange than for ir + echange, and (c) most importantly the two thresholds \D°D*° > and 
\D + D*~ > are spit by 8.1 MeV. We still assume that the couplings in the pion exchange are related by exact isospin 
which gave the matrix (5). Clearly one expects that there will be only one bound state (i.e. if there is a bound state 
near \D a D* a > there should be no state near \D + D*~ > nor ID^D* >). The perturbation from isospin breaking 
can only increase the binding of the lower state and make the orthogonal state feel even more repulsion. 

One then has a four coupled channels, two spin orbitals (S and D-wave) coupled with two charge modes \D°D*° > 
and the \D + D*~ >. The computed wave functions are shown in Fig. 3 when one assumes a binding energy of half 
an MeV with respect to the lower threshold \D°D*° >. 
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FIG. 3. The 4 different components of the wavefunction for a 1 ++ DD* deuson when isospin is brokren by the isospin mass 
splittings of the D, D* and n. The binding energy with respect to the D°D*° threshold is assumed to be 0.5 MeV. Note 
the strong dominance of the \D°D*° > (full lines) compared to the \D + D*~ > (dashed) in the S-wave. (The regularization 
parameter A is here 1.40 GeV). 

One notices that the \D a D*° > component in the wave function strongly dominates over the \D + D*~ > in the 
S- waves for larger distances as could be expected. Therefore in terms of isospin states there is a substantial 1=1 
component in the state. As a consequence of this the decay of the X(3872) to p J/t/j — ► ir + ir~ J/ip will not be 
forbidden. There are indications of this in the experiment [1]. On the other hand the decay chain aJ/tp — ► tt + it~ J/i/j, 
where a is any isoscalar object, is forbidden by C-parity and spin-parity for a J PC = ^ or 1 ++ deuson. Thus one 
should expect some p°J/ip — > 7r + 7r _ J/ip, but no aJ/ip and in particular no J/ipir a ir°. An experimental observation of 
the J/ifjir°ir a would exclude the X(3872) to be a deuson. This is unfortunately a difficult mode to see experimentally 
and it seems easier to with more statistics confirm that the tt + tt~ in ir + n~ J/ip comes only from p. 

The second case when pion exchange is attractive and strong is for a P-wave (near) isoscalar DD* system with 
J PC = The important thing here to notice is that the large tensor term adds directly to the central term and 
not as in eq.(4) through off-diagonal terms and iteration. The large 1/r 3 term can then in principle overcome the 
P-wave angular momentum barrier! The potential is 
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V -+ = lV [C{r) + 2T(r)] = -ji±V 
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(6) 
(7) 



which after a regularization can bind a pseudo-scalar DD* deuson in spite of the P-wave. The existence of this state 
is however more model dependent than in the previous case since it is much more sensitive to the regularization of 
the tensor force, and to other short range contributions. 

Thus although the quantum numbers J PC — 1 ++ seems favoured J PC = 0"+ is still not excluded that the X(3872) 
is a pseudoscalar DD* deuson. The fact that in that case the central and the large tensor potential adds directly 
to the same spin orbital is similar to a two nucleon system with 3 Pq quantum numbers, but with the important 
difference that one now has attraction, not repulsion from pion exchange. In Fig. 4 (assuming exact isospin) and Fig. 
5 (including isospin breaking) we show the wave function for such a possible pseudoscalar deuson when the binding 
energy is again set to be 0.5 MeV below the lower threshold. 
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FIG. 4. The DD* pseudoscalar deuson wave function in the exact isospin limit, when the binding energy is assumed to be 
0.5 MeV. (The regularization parameter A is here 1.46 GeV). 
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FIG. 5. The \D°D*° > and \D + D*~ > pseudoscalar deuson wave functions with isospin breaking and when the binding 
energy below the D°D*° threshols is 0.5 MeV. (The regularization parameter A is here 1.53 GeV.) 

Concluding remarks 

We conclude with a few general comments. The heavier the constituents are the stronger is the binding, since the 
kinetic repulsion becomes smaller and is more easily overcome by the attraction from the potential term. Thus as 
seen from table 1 the DD* and D*D* systems are barely bound but for BB* and B*B* the binding energy is w 50 
MeV [21]. 
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An uncertainty in the original calculation [21] was the D* coupling to Dir, which was modelled from the NNir 
coupling. We predicted a D* + — > D°ir + width of 63.3 keV in excellent agreement with the recent measurement of 
65 ±3 keV [25]. This increases the reliability of that calculation, and we have here kept the same overall normalization 
V = 1.3 McV. 

For flavour exotic two-meson systems (/ = 2, double strange, charm or bottom), such as DD* or B*B* } pion 
exchange is always either too weakly attractive or repulsive. Calculations do not support such bound states to exist 
from pion exchange alone, and shorter range forces are expected to be repulsive. Should flavour exotic BB* or B*B* 
states exist, however (See [23]), they would be quite extremely narrow since they would be stable against strong 
decays. 

In this calculation we have for simplicity not included the threshold effects from u)J/ip, p°J/ip nor three body 
channels DDtt etc. into the coupled channel framework. In particular u)J/ip (which is just above the X(3872)) can 
be important. Swanson [12] finds a rather large coupling of the X(3872) to coJ/ip when including a short range 
quark exchange potential. Then a substantial coJ/ip component can be generated into the wave function. This then 
contributes to the binding of X(3872) through higher order terms like X — > u>J/tp — ► X. This mechanism is similar 
to the fact that the deuteron binding increases because of the freedom for the 3 5i state to turn into a 3 Di and back. 
It would be important to have a measurement of the ratio of branching ratios, B(X — ► D a D°Tr a )/B(X — ► J/^tttt), 
which should not be too small for X as a deuson. Belle [26] has measured an upper limit for the product B(B + — > 
XK + ) x B(X — » L»°Z) 7r ) < 6 x 10~ 5 , but this is not yet significant because of the small B(B + — > XK + ) compared 
to B{B+ -» 7/>(3770)iV+). 

Another recently seen peak where pion exchange should be important is the observation of BES [27] in J/ip — ► jpp 
near the pp threshold. Pion exchange should be the most attractive for TV TV quantum numbers when 1=0, L=0, J=0 
or for a pseudoscalar 1 S state. These are also the likely quantum numbers for the peak which they have observed 
just above the pp threshold. However, if real, it need not necessarily be due to a resonance below theshold, but could 
be due to a strong final state interaction involving pion exchange between the produced proton-antiproton system 
[28] . Similar pp peaks but with low statistics have also been seen in B ± -» pp^ [29] and B° ppD^° [30] 

A more detailed understanding with further experimental information on the Belle charmonium state is important. 
If our arguments are supported by data it could open up a completely new spectroscopy for heavy mesons. It might 
also throw some new light on many problematic light resonances in particular the /i(1420), 7/(1410) and 7/(1480) 
resonances just above the KK* threshold, where pion exchange in the final state interaction can play an important 
role, although it is not expected to be strong enough to alone bind those states. Possibly even the glueball candidate 
/o(1500) could have a large (pp — u)cu)/\/2 component due to the strongly attractive pion exchange for those quantum 
numbers [20]. 

Acknowledgement I thank Steve Olsen and San Fu Tuan for useful comments 
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